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Abstract 

Magnetization reversal in magnetic nanostructures is investigated numerically over 
time-scales ranging from fast switching processes on a picosecond scale to thermally 
activated reversal on a microsecond time-scale. A simulation of the stochastic 
Landau-Lifshitz equation of motion is used as well as a time quantified Monte 
Carlo method for the simulation of classical spin systems modeling magnetic Co 
nanoparticles. For field pulses larger than the Stoner-Wohlfarth limit spin pre- 
cession effects govern the reversal behavior of the particle while for lower fields a 
magnetization reversal is only possible when it is assisted by thermal fluctuations. 
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1 Introduction 



The miniaturization of magnetic structures plays an important role for funda- 
mental research as well as for technical applications. This leads to an incremen- 
tal interest in the understanding of the behavior of small magnetic particles and 
structures down to the nanometer scale. But with decreasing size of the magnetic 
system thermal activation becomes relevant. Hence, much effort is focused now on 
the understanding of magnetization dynamics at finite temperatures. 

In the following, we briefly describe numerical techniques for the study of mag- 
netization dynamics in nanostructures, modeled as classical spin systems, where 
a finite temperature is taken into account. First, we will focus on the underlying 
model and the two basic methods, Langevin dynamics and Monte Carlo simu- 
lation . Then we study the magnetization reversal in Co nanoparticles. Starting 
with the deterministic spin dynamics on short time scales which plays a crucial role 
in high-speed data storage we go on to the probabilistic long-time behavior 

where a thermally assisted reversal can occur even for magnetic fields below the 
coercive field HTHTlll. 



2 Classical Spin Model 

The micromagnetic properties of a nanoparticle can be described using a model 
of classical magnetic moments which are localized on a given lattice. Such a spin 
model can be motivated following different lines: on the one hand it is the classical 



limit of a quantum mechanical, localized spin model — the Heisenberg model [|12 
On the other hand, a classical spin model can also be interpreted as the discretized 
version of a micromagnetic continuum model |13[, where the charge distribution 



for a single cell of the discretized lattice is approximated by a point dipole |T^,p!4 

The interpretation as an atomic model restricts the use of computer simulations 
to the investigation of rather small systems of only a few million atoms - corre- 
sponding to particle sizes of only a few nanometers. On the other hand, within a 
continuum model, the space might be discretized on a much larger length scale, 
as compared to an atomic distance. However, in continuum theory usually a con- 
stant absolute value of the magnetization vector is assumed, an assumption which 
fails for higher temperatures since the space averaged magnetization breaks down 
when approaching the critical temperature. Hence, one expects correct thermal 
properties only in the limit of small cell sizes of the order of atomic distances. 
In the following, let us consider a classical three dimensional Heisenberg Hamil- 
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tonian for localized spins, 



n = -j^s,-s,-/i,B-^s,-4E('5r)' 

(1) 

where the Sj = n^/ fis are three dimensional magnetic moments of unit length rep- 
resenting atomic magnetic moments. The first sum is the ferromagnetic exchange 
of the moments with the coupling constant J. The second sum is the coupling of 
the magnetic moments to an external magnetic field B, the third sum represents a 
uniaxial anisotropy, here, for dz > favoring the z axis as easy axis of the system, 
and the last sum is the dipolar interaction where w = HqhH {A-na?) describes the 
strength of the dipole-dipole interaction. The are unit vectors pointing from 
lattice site i to j and is the distance between these lattice sites in units of a. The 
dipole-dipole interaction can be computed efficiently with the help of fast Fourier 
transformation methods |TB|,|TB|. One should however note that in a Monte Carlo 
simulation with a single-spin flip algorithm the FFT method is an approximation 
the implementation of which was described in details before ITT 



The transformation of the above introduced atomic parameters to the material 
parameters which are usually used in a continuum model is given by J = 2aAx 
where is the exchange energy, dz = Ka^ where K is the anisotropy energy 
density, and /x^ = MgO? where Mg is the spontaneous magnetization. 



3 Landau-Lifshitz Equation and Langevin Dynamics 

In the short time limit spin precession is important which can be taken care 
of by studying the corresponding equation of motion. The basic numerical ap- 
proach which includes thermal activation, is the direct numerical integration of 
the Langevin equation of the problem. In order to obtain thermal averages one 
has to calculate many of these trajectories starting with the same initial condi- 
tions taking an average over these trajectories for the quantities of interest. This 
method is referred to as the Langevin dynamics formalism 0. 

The underlying equation of motion for a magnetic system is the Landau- 
Lifshitz-Gilbert (LLG) equation. 



dt {l + a^)^ls 



H,(t) + a(S, X H,(t)) 



with the gyromagnetic ratio 7 = L76 x lO^^(Ts) ^, the dimensionless damping 
constant a, and the internal field Hj(t) = Ciit) ~ dH/dSi. Langevin dynamics is 
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introduced here in form of the noise Ci(^) which represents thermal fluctuations, 
with {Ciit)) = and (C(^)Cj(^')) = '^^ij^-ne^it — t')akBTfis/l where i,j denote 
once again lattice sites and r], 9 Cartesian components. 

The LLG equation with Langevin dynamics is a stochastic differential equation 
with multiplicative noise. For this kind of differential equation a problem arises 
which is called the Ito-Stratonovich dilemma ITT]. As a consequence, different time 



discretization schemes may converge to different results with decreasing time step. 
As was pointed out in |18[ the multiplicative noise in the Langevin equation above 



has to be treated by means of the Stratonovich interpretation. An appropriate 
discretization scheme leading to a Stratonovich interpretation is the Heun method 
Ijl, 17, 1^ which is used in the following. 



4 Monte Carlo Methods 

In the long time limit only spin relaxation and thermal fluctuations are relevant 
which can be studied very conveniently using Monte Carlo methods with quantified 
time step 0. Within a Monte Carlo approach trajectories in phase space 
are calculated following a master equation ||2y] for the time development of the 
probability distribution Pgif) in phase space, 

dPs 

= Y^iPs'Ws'-^s - PgWs^s')- (3) 

Here, s and s' denote different states of the system and Ws'^s is the transition rate 
for a change from a state s' to a state s. These rates have to fulfill the condition pOl 



exp 



E{S) - E{S') 



(4) 



The master equation describes exclusively the coupling of the system to the heat 



bath |]20|. Hence, only the irreversible part of the dynamics of the system is 
considered including only the relaxation and the fluctuations. A Monte Carlo 
simulation does not include the energy conserving part of the equation of motion. 
Hence, no precession of magnetic moments will be found. 

Monte Carlo approaches in general have no physical time associated with each 
step of the algorithm. However, recently a time quantified Monte Carlo method was 



proposed in and later successfully applied to different model systems [[10| , pT| ,pi 
Here, the interpretation of a Monte Carlo step as a realistic time interval At was 
achieved by a comparison of one step of the Monte Carlo process with a time 
interval of the LLG equation in the high damping limit. We will use this algorithm 
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in the following. The trial step of this algorithm is a random movement of the 
magnetic moment within a cone with a given size r with 



2 20A;BTa7 



(1 + a^)Hs 

In order to achieve this efficiently one constructs a random vector with constant 
probability distribution within a sphere of radius r. This random vector is added 
to the initial moment and subsequently the resulting vector is normalized 0. 

Using this algorithm one Monte Carlo step represents a given time interval At 
of the LLG equation in the high damping limit as long as At is chosen appropriately 
(for details see [0]). 



5 Precessional Reversal in Co Nanoparticles 

In the following we consider Co nanoparticles, where the material parameters 
are = 1.3 ■ 10"" 3/m, K = 6.8 ■ 10^ J/m^ and M, = 1.43 ■ 10*^ A/m. For 
simplicity we simulate a simple cubic lattice with atomic distance a = 0.25 nm. 
Our simulation starts with a spin configuration where all magnetic moments point 
into the z direction, aligned with the easy axis, and with the z component of the 
external magnetic field antiparallel to the magnetization so that the system is in 
an unstable, or at least metastable state. 

In sufficiently small particles the magnetic moments rotate coherently during 
the magnetization reversal. A quantitative description of coherent rotation in 



ellipsoidal single domain particles was developed by Stoner and Wohlfarth [|2 
Depending on the angle between the applied field B and the z (easy) axis of 
the system, the coercive field Be varies following the so-called Stoner- Wohlfarth 
asteroid [^]. Under an angle of 45° it is Be = dcsV/fig where V is the volume of 
the particle and des is an effective anisotropy constant. In our case it is 5c ~ 0.7 
T. 

First, we are interested in fast switching processes where the applied field is 
higher than the coercive field and the reversal is dominated by spin precession. We 
simulate the system with Langevin dynamics as described before. Figure |l| shows 
the time dependence of the magnetization of our Co particle in the low damping 
limit {a = 0.1) for a very small ellipsoidal shaped particle of length L = 4nm and 
diameter L = 2nm. 

The z and the x components of the magnetization are shown as well as its 
absolute value. The magnetic field B is set under an angle of 45° to the z axis 
within the yz plane so that the response of the system to the external field sets 
in directly. The wavering magnetization of the system clearly follows from the 
precession of the spins. 
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Figure 1: Reduced magnetization vs. time for a Co nanoparticle. The data are 
from Langevin dynamics simulations. B = 1.13 T and T = 16 K. 

Note, that the precession time of our system is not simply given by the pre- 
cession time of a single spin in an external field (xp = 27r(l + a'^)/'~fB ^ 32 ps in 
our case). Instead, the whole internal field is relevant for the spin precession, also 
the contributions from the dipolar field, the exchange and the anisotropy. This 
internal field is not constant in time and in general it is non-homogenous within 
the system. However, in the very small particle which we consider here, the in- 
ternal field is sufficiently homogenous so that the magnetization moves coherently. 
Hence, the absolute value of the magnetization remains constant in time as is 
shown in the figure. Note also, that even after the new stable state is reached 
the magnetization still keeps on oscillating around its equilibrium value, driven by 
thermal fluctuations. 

6 Thermally Activated Reversal 

Let us now turn to the case B < Be, where the reversal process can only occur 
when it is thermally activated. Since we are now interested in the long time and 
high damping limit {a = 4) where the behavior of the particle is governed by 
thermal fluctuations we can use Monte Carlo methods. Figure ^ shows the typical 
time dependence of the magnetization of the same Co particle as before. The field 
B is set antiparallel to the initial state here, so that the zero-temperature coercive 
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Figure 2: Reduced magnetization vs. time for the same particle as in Fig. [1|. The 
data are from Monte Carlo simulations. 5 = I.IT and T = 5 K. 



field is given by = 2dcsV/ /ig which in our case is ~ 1.4 T. 

As one can see, the z component of the magnetization remains nearly constant 
for a time period which is rather long as compared to the previous simulation. 
Then, suddenly, the magnetization drops and its sign changes. As before the ab- 
solute value of the magnetization is constant in time and the reversal mechanism 
is mainly a coherent rotation. The value of the switching time is approximately 
6.6 ns in our simulation. However, this thermally activated switching is not a 
deterministic process as it was the short time dynamics studied before, where the 
switching followed mainly from the deterministic part of the equation of motion. 
Instead, the thermal activation process here is a probabilistic event. The proba- 
bility distribution P{ts) for switching events taking place after a time tg follows 
an exponential law 0, 

P(t,) ~exp(-ts/r), (6) 
in the limit of large time-scales. Here, r is a characteristic time scale. 



r = Toexp (AE/kBT), 



(7) 



where tq is a prefactor and AE an energy barrier which both are related to a certain 
reversal mechanism (see e. g. [0,^^ for analytically determined prefactors and 
energy barriers in different systems and [^^,[1^, ^ for numerical work on this 
subject). In general the prefactor may depend on the system parameters, the 
temperature, the applied magnetic field and the damping constant. 
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Figure 3: Characteristic time r vs. 1/T for the same Co particle as in the figures 
before. The slope of the solid line represents the energy barrier AE. B = 1.1 T. 

For the case of a Stoner-Wohlfarth particle with the applied field parallel to 
the easy axis the energy barrier has the form 

AE = d,f,v(l-^y. (8) 

-De 

This energy barrier as well as the prefactor have been calculated by Brown [0 
under the assumption that all magnetic moments are parallel, so that the system 
behaves like one single magnetic moment. 

For a further analysis we extract the energy barrier which governs the reversal 
process from our numerical data. Figure ^ shows the temperature dependence of 
the characteristic time, i. e. the mean switching time, obtained from our simula- 
tions. The slope of the solid line corresponds to the theoretical value of the energy 
barrier for a reversal by coherent rotation obtained from Eq. |^. Obviously, it is in 
very good agreement with our numerical data for low enough temperatures. 
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